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ABSTRACT 

We present analysis of the UV-spectrum of the \ow-z AGN IRAS-F22456- 
5125 obtained with the Cosmic Origins Spectrograph on board the Hubble Space 
Telescope. The spectrum reveals six main kinematic components, spanning a 
range of velocities of up to 800 km s" 1 , which for the first time are observed in 
troughs associated with C II, C IV, N v, Si II, Si ill, Si IV and S IV. We also obtain 
data on the O VI troughs, which we compare to those available from an earlier 
FUSE epoch. Column densities measured from these ions allow us to derive 
a well-constrained photoionization solution for each outflow component. Two 
of these kinematic components show troughs associated with transitions from 
excited states of Sill and Cll. The number density inferred from these troughs, 
in combination with the deduced ioinization parameter, allows us to determine 
the distance to these outflow components from the central source. We find these 
components to be at a distance of ~ 10 kpc. The distances and the number 
densities derived are consistent with the outflow being part of a galactic wind. 

Subject headings: galaxies: quasars — galaxies: individual (IRAS F22456-5125) - 
line: formation — quasars: absorption lines 
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INTRODUCTION 



Mass outflows are detected in the UV spectra of more than 50% of low redshift act ive 
galact ic nuclei (AGN) ma i nly Seyfert galaxies, e.g . 



( 120021 ) 



Dunn et al 



(120071 ) 



Crenshaw et al 



(1999) 



Kriss et al. 



Ganguly fc Brothertonl (120081 ) . These outflows are observed as 



narrow absorption lines (a few hundred km s 1 in width) blueshifted with respect to the 
AGN systemic redshift. 

In this paper, we determine the ionization equilibrium, distance, mass flow rate, and 
kinetic luminosity of the UV outflow observ ed in the luminous Seyfert 1 galaxy IRAS 



F22456-5125 (z = 0.1016. iDunn et al 



2010l ). The bolometric luminosity of this object, 



^bol 



10 45 6 ergs s 1 (see Section @J, places it at 



;he Seyfert /quasa r border defined to 



be 10 I/©, where L is the luminosity of the sun ( iSoifer et al. 



19871 ). Several absorption 



systems are resolved in the UV spectrum in five main kinematic components ranging in 
velocities from —20 km s _1 to —820 km s _1 . A detailed analysis of the physical properties of 
the UV absorber determined from Far Ultraviole t Spectroscopic Explorer (FUSE) archival 
spectra has been published by lDunn et al.l (120101 ) . These authors report a lower limit on the 
distance R of the absorbing material from the central source of 20 kpc using photoionization 
timescale arguments. 

In June 2010 we observed IRAS F22456-5125 with the Cosmic Origins Spectrograph 
(COS) on board the Hubble Space Telescope (HST) as part of our program aiming at 
determining the cosmological impact of AGN outflows (PI: Arav, PID: 11686). The high 
signal to noise spectrum obtained reveals the presence of absorption troughs associated with 
high ionization species (Civ, Nv, O VI, Si IV and Siv) as well as lower ones (Sill, Si ill, Cll) 
thus increasi ng the numb e r of c onstraints on the photoionization analysis of the absorber 



compared to 



Dunn et al. 



(120101 ). We also identify absorption troughs corresponding to 
excited states of Si II and C II associated with two kinematic components of the UV outflow. 
The population of the excited state relative to the resonance counter part provides an indirec t 



measurement of the number density of the gas producing the lines (lOsterbrock fc Ferland 
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20061 ). These number densities allow us to determine reliable distances to these two 
components and hence derive their mass flow rates and kinetic luminosities. 

The plan of the paper is as follows: in § [2] we present the COS observations of IRAS 
F22456-5125 as well as the reduction of the data and identification of the spectral features 
within the COS range. In §[3] we detail the computation of the column densities associated 
with every species. We present the photoionization analysis of the outflow components 
in § H] and report the derived distance, mass flow rate, and kinetic luminosity in § [5j We 



conclude the paper by a discussion of our results in 



series and the reader will be referred to 



61 Th is paper is the second of a 



Edmonds et al. 



(12011 



hereafter Paper I) for further 



details on the techniques used throughout the paper. 



2. HST/COS observations and data reduction 



We observed IRAS F22456-5125 using the COS instrument (lOsterman et al.ll2010l ) on 
board the HST in June 2010 using both medium resolution (AA/A ~ 18, 000) Far Ultraviolet 
gratings G130M and G160M. Sub-exposures of the target were obtained for each grating 
through the Primary Science Aperture (PSA) using different central wavelength settings in 
order to minimize the impact of the instrumental features as well as to fill the gap between 
detector segments providing a continuous coverage over the spectral range between roughly 
1135-1795 A. We obtained a total integration time of 15,056 s and 11,889 s for the G130M 
and G160M gratings, respectively. 



>fl pil 



The datasets processed through the standard CALCOqj pipeline were retrieved 
from the MAST archive. They were then flat-fielded and combined together using the 



COADDJQDB IDL pipeline developed by the COS GTO team (see 



Danforth et al 



2010 



for details). The reduced spectrum with its original ~ 2 km s 1 oversampling has an overall 



1 Details on CALCOS can be found in the COS Data Handbook. 

2 The routine can be found at http:/ /casa. colorado.edu/ ~danforth/science/cos/costools.html 
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signal to noise ratio > 15 per pixel in most of the continuum region. Typical errors in the 
wavelength calibration are less than 15 km s _1 . In Figured], we show the majority of the 
spectrum on which we identified major intrinsic absorption features associated with the 
outflow. The COS FUV spectrum of IRAS F22456-5125 is presented in greater detail along 
with the identification of most absorption features (interstellar, intergalactic, and intrinsic 
lines) in the on-line version of Figure [TJ 



2.1. Identification of spectral features 



Using archived FUSE spectra iDunn et al.l (120071 . |2010| ) reported the first detection 



of five distinct kinematic components with centroid velocities v\ = —800 km s 1 , 
v 2 = -610 km s _1 , v 3 = -440 km s"\ v 4 = -320 km s _1 , v 5 = -130 km s _1 , and FWHM 
e [50, 200] km s" 1 associated with an intrinsic UV outflow in IRAS F22456-5125. These 
components, spanning a total velocity range of 800 km s -1 , were detected in O VI, Cm and 
in several lines of the Lyman series (Ly/3 to Lyr?). Using the kinematic pattern reported 



by 



Dunn et al. 



(120101 ) as a template we identify absorption features in our COS spectrum 
related to both low ionization (Cll, Sill, Sim) and high ionization species (Si IV, SlV, Civ, 
N v, O vi) as well as in the Lya transition. Absorption troughs from the metastable level 
Cll* A1335.704 are detected in components 2 and 3, and troughs from metastable Sill* 
A1264.738 and A1194.500 are detected in component 2. 

While the absorption troughs associated with the higher ionization lines generally 
exhibit broader profiles, we observe a 1:1 kinematic correspondence between the core of 
these components and the narrower features associated with the lower ionization species of 
the outflow. Given the significantly broad range of velocities covered by the components and 
their net kinematic separations, such a match is not likely to occur by chance. This argues 
in favor of a scenario where the troughs of the different ionic species detected in a given 
kinematic component are generated in the same region. This observation is strengthened 
by the fact that most of the troughs have a line profile similar to that of the non-blended 
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Fig. 1.— The full FUV spectrum of IRAS F22456-5125 obtained by COS. The major ab- 
sorption troughs related to the intrinsic absorber are labeled. A full identification of all the 
absorption features is presented in the online version. The green line represent our fit to the 
non absorbed emission model (see Section l2~3i) . 
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Nv A1238.820 line when properly scaled. 

The high S/N of our COS observations (S/N > 40 per resolution element on most of 
the spectral coverage) reveals the presence of kinematic substructures in several component s 



of the outflow compared to the lower S/N FUSE observations (S/N ~ 7. lDunn et al.ll2010l ). 
Nevertheless given the self-blending of these features in the strongest lines (e.g. O vi) and 
the absence of apparent change betw een the FUSE and COS observations, we will use the 



labeling of components as defined in 



Dunn et al. 



(120101 ) . We will however separate their 



trough 5 into low and high velocity components given the apparent difference in ionization 
suggested by the presence of a stronger Si ill in subcomponent 5A {v§a — ~ 40 km s -1 ) 
than in 5B (v§b — —130 km s _1 ) relative to the higher ionization lines (Civ, Nv, O VI, see 
Figure [3]). Most of our analysis in this paper concentrates on components 2 and 3 of the 
outflow, for which absorption features associated with an excited state have been detected. 



2.2. Deconvolution of the COS spectrum 

Detailed analysis of the on-orbit COS Line Spread Function (LSF) revealed the 
presence of broadened wings that scatter a significant part of the continuum flux inside the 



absorption troughs (see 



Kriss et al 



2011 



for details). This continuum leaking is particularly 
strong for narrow absorption troughs (FWHM ~ 50 km s -1 ) in which this effect may 
significantly affect the estimation of the true column density by artificially increasing the 
residual intensity observed inside the troughs. 

Given the overall good signal to noise ratio of our data, we can correct the effect 
of the poor LSF by deconvolving the COS spectrum. Adopting the procedure described 



in 



Kriss et al.l ( 1201 ll ). we deconvolve the spectrum obtained for each grating in 50 



A intervals using the wavelength dependent LSFs and an IDL implementation of the stsdas 
Richardson-Lucy "lucy" algorithm (G. Schneider & B. Stobie, private communication, 
2011). The main effect of the deconvolution is illustrated in Figure [2] in which we clearly 
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see that the deconvolved spectrum shows significantly deepened intrinsic Lya absorption 
troughs as well as produces a square, black bottom for the saturated interstellar line C II 
A1334.532. 

However, the main drawback of the deconvolution algorithms commonly used, such as 
the Richardson-Lucy (RL) algorithm, is a significant increase of the noise in the deconvolved 
spectrum due to the fact that these techniques try to perform a total deconvolution of 
the data, i.e. in which the LSF of the decon yolved spectrum is a Dirac delta function, 



violating the Shannon sampling theorem (see 



Magain et al. 



19981 for a thorough discussion 



of these issues). In order to decrease these effects, we modified the RL algorithm by forcing 
the deconvolved spectrum to have a LSF satisfying the sampling theorem. We choose the 
deconvolved LSF to be a Gaussian with a 2 pixel FWHM (~ 5 km s -1 given the COS 
detector sampling). This operation prevents the appearance of strong unwanted oscillations 
since we force the maximum resolution that can be achieved in the deconvolved data to 
agree with the sampling theorem. The deconvolved spectrum produced by this modified 
Richardson Lucy algorithm is similar to the one obtained by the traditional RL algorithm 
(see Figure [2]), the main difference being the significant decrease of the high frequencies and 
high amplitude features artificially generated by RL deconvolution with a high number of 
iterations. In our analysis, we will derive the column density for each ionic species using 
the spectrum deconvolved with the modified RL algorithm, allowing us to derive more 
accurate column densities associated with the narrow absorption components observed in 
IRAS F22456-5125. 



2.3. Unabsorbed emission model 

The unabsorbed emission model ^(A) of IRAS F22456-5125 is constructed in a 
similar manner to the one described in detail for IRAS-F04250-5718 in Paper I, in which 
we consider three main sources of emission : a continuum, a broad emission line (BEL) 
component and a narrow emission line (NEL) component. Adopting a single power law 
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Fig. 2. — Illustration of the necessity of using a deconvolution algorithm when dealing with 
COS data (see text for details). Troughs associated with the intrinsic absorber significantly 
deepen while the saturated interstellar C II line exhibits the expected squared black bottom 
profile. The main difference between the deconvolved spectrum using the RL method and 
the modified RL algorithm respecting the sampling theorem is the significant reduction of 
oscillations due to the total deconvolution process performed in RL, even more so when 
considering a high number of iterations. 
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F(X) = Fii 50 (A/1150) a to describe the deredenned (E(B-V)=0.01035. ISchlegel et al.lll998f l 
continuum emission, we obtain a reduced Xred = 1-413 over emission/absorption line free 
regions of the rest wavelength spectrum ([1115, 1130] A, [1340,1360] A, and [1455,1475] A) 
with a = -1.473 ± 0.068 and Fn 50 = 2.130 10" 14 ± 0.0033 10~ 14 erg cm^A-V 1 . 

Prominent BEL features observed in the spectrum (Lya, C IV, O vi) are fit using two 
broad gaussians of FWHM ~ 9000 and 2000 km s . The NEL component of each line of 
a doublet is fit by a single narrower gaussian (FWHM ~ 600 km s _1 ) centered around the 
rest wavelength of each line, with the separation of the two gaussians fixed to the velocity 
difference between the doublet lines. The NEL of the strong Lya line is best fit by two 
gaussians of FWHM ~ 1200 and 400 km s _1 . The remaining weaker emission features in 
the spectrum (SilV+OlV, Cn, Nv, Ol etc.) are modeled by a smooth cubic spline fit. A 
normalized spectrum is then obtained by dividing the data with the emission model. We 
present our best fit to the unabsorbed spectrum of IRAS F22456-5125 in Figured] 



3. Column Density Determination 



3.1. Methodology 



The column density associated with a given ionic species detected in the outflow is 
determined by modeling the residual intensity in the normalized data of the absorption 
troughs. Assuming a single homogeneous emission source F (v) and a one-dimensional 
spatial distribution of optical depth ac ross the emission source Ti(x, v), we can express the 



intensity Fi(v) observed for a line i as (lArav et al.ll2005l ): 



FAv) 



F (v) I < 
lo 



i(x,v) d 



X 



where v is the radial velocity of the outflow, and the spatial extension of the emission source 
is normalized to 1. Once the optical depth solution Ti(x, v) is derived at a given radial 
velocity, we link the observed residual intensity Ii(v) = Fi(v)/Fo(v) to the ionic column 
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density using the relation : 

N ion (v) 



3.8 x 10 



11 



< Ti(v) > (cm 2 km 1 s) 



(2) 



where /j, A i; and < Ti(v) > are the oscillator strength, the rest wavelength, and the average 
optical depth across the emission source of line i (see Paper I), respectively. 

We consider here the three absorber models (i.e. optical depth distributions) discussed 
in Paper I; the Apparent Optical Depth (AOD), Partial Covering (PC), and Power Law 
(PL) models. We use these three models in order to account for possible inhomogeneities 
in the absorber (see Sect. |6]), which cause the apparent strength ratio R a = Ti/rj of two 
lines i, ? from a given ion to deviate from the expected la boratory ratio Ri = \ifi/\jfj 



(e.g. 



Wampler et al 



1995 



Hamann 



1997 



Arav et al 



1999). Wherever possible we derive 



these three optical depth solutions for ions with multiple transitions. However, as mentioned 
in Paper I, we consider the results obtained with the PL model performed on doublets with 
caution given its increased sensitivi ty to the S/N, wh ich can lead to severe overestimation 



of the underlying column density ( lArav et al 



20051 ). For singlet lines we will generally 



only derive a lower limit on the column density using the AOD method. This lower limit 
will be considered a measurement in cases where the singlet line is associated with a 
kinematic component for which other multiplets do not show signs of saturation. In the 
following subsections, we use the term (non-black) saturation to qualify absorption troughs 
of doublets in which R a = Ti/rj < 0.75 R[, where Tj and Tj are the apparent optical depth 
of the strongest and the weakest component of the doublet, respectively. 



3.2. Column density measurements 

Computed ionic column densities are determined using the deconvolved line profiles 
presented in Figure Eland the ionic transition properties reported in Tabled! The computed 
column densities are reported in Table [2] for the three absorber models when possible. The 
adopted values shown in the last column of Table |2] are the ones used in the photoionization 
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analysis. When available, we choose to use the value reported in the PC column as the 
measurement and use the PL measurement and error as the upper error in order to account 
for the possible inhomogeneities in the absorbing material distribution. If only the AOD 
determination is available we will consider the reported value minus the error as a lower 
limit unless we have evidence suggesting a high covering. 



3.2.1. Hi 



The spectral coverage of the COS G130M and G160M gratings only allows us to 
cover the Lya line that shows a deep and smooth profile in which the different kinematic 
components blend. The absence of higher-order Lyman series lines restrict us to put a lower 
limit on the H I column density by applying the AOD method to the Lyct profile. 

A better constraint on the H I col umn density is det ermined by using higher-order 



Ly-series lines from earlier FUSE data 



Dunn et al. 



( 120101 ) . In Figure H] we compare the 



June 2010 HST/COS and 2004 FUSE spectra of IRAS F22456-5125 in the overlapping 
region, essentially showing the rest frame O VI region of the spectrum. While we observe a 
net increase in the continuum flux (a factor of ~ 1.17 between the higher S/N 2002 FUSE 
observation and the 2010 COS observation) and in the broad emission line flux between the 
two epochs, the overall shape of the absorption troughs and continuum remains unchanged. 
Careful examination of the normalized O VI COS and FUSE absorption line profiles reveals 
troughs that are consistent, given the limited S/N of the FUSE observations, with no 
variations between the epochs in any of the kinematic components. Therefore we use the 
Hi column densities determined in the FUSE spectrum i n our analys i s. Th e Hi column 



Dunn et al. 



(120101 ) and consist 



density estimates reported in Table H] are extracted from 
of partial covering solutions derived on higher-order Ly-series lines. We however note 
report the detection of the LylO line associated with kinematic c omponent T2 reve aling an 



underestimation of the Hi column density by a factor of ~ 2 in 



Dunn et al. 



(120101 ) . The 



AOD solutions we report are computed on the Lya line profile present in the COS spectrum 
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Fig. 3. — Normalized absorption line profile of the metal lines associated with the outflow in 
IRAS F22456-5125. The line profiles have been deconvolved using the modified RL algorithm 
described in Section 12.21 and rebinned to a common ~ 5 km s -1 dispersion velocity scale. 
For doublets, we overplot the expected residual intensity in the strongest component based 
on the residual flux observed in the weakest component assuming an AOD absorber model. 
For ClV we only plot that quantity in regions free of self blending (mainly T3, see text). 
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Table 1. Atomic Data for the Observed Transitions 



Ion 


Z7 a 




^low 


ft 




cm -1 


A 






Hi 


0.00 


1215.670 


2 


0.4164 


Cii 


0.00 


1334.532 


2 


0.1290 


Cll* 


63.42 


1335.704° 


4 


0.1277 


Civ 


0.00 


1548.202 


2 


0.1900 


Civ 


0.00 


1550.774 


2 


0.0952 


Nv 


0.00 


1238.821 


2 


0.1560 


Nv 


0.00 


1242.804 


2 


0.0780 


Ovi 


0.00 


1031.912 


2 


0.1330 


Ovi 


0.00 


1037.613 


2 


0.0660 


Sin 


0.00 


1190.416 


2 


0.2770 


Sin 


0.00 


1193.280 


2 


0.5750 


Si II* 


OCT O A 

287.24 


1194.500 


4 


A TO TA 

0.7370 


Sin 


0.00 


1260.422 


2 


1.2200 


Sin* 


287.24 


1264.730 


4 


1.0900 


Sin 


0.00 


1304.370 


2 


0.0928 


Sin 


0.00 


1526.720 


2 


0.1330 


Sim 


0.00 


1206.500 


1 


1.6700 


Si iv 


0.00 


1393.760 


2 


0.5130 


Si iv 


0.00 


1402.770 


2 


0.2550 


Siv 


0.00 


1062.656 


2 


0.0500 



a - Lower level energy, b - Wavelength of the transition, c - Statistical weight, d - Oscillator strength. We 
use the oscillator strengths from the National Insti tute of Standards and Technology (NIST) database, 



except for S IV for which we use the value reported in 



Hibbert et al 



(|2002l ). e - Blend of two transitions, we 



report the sum of the oscillator strength and the weighted average of Aj 
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Table 2. Computed column densities 



Trough Vf Ion AOD a PC" PL a Adopted*' 

km s -1 10 12 cm~ 2 10 12 cm~ 2 10 12 cm~ 2 10 12 cm~ 2 



Tl 


-800 


II I 

Civ 


73 9+ 04 
20.9+11 


< 900 b 




6 [73.5, 900] 
20.9+. 11 , 






N v 






20.2+*;<j 








Ovi 


474^ 






745+i 3 / 






Sim 


< °-n±° °oi 






< 0.36 


T2 


-610 


H i 


436ti 6 


4*>°tSSS b 




£ [9400, 15800°] 






Cn 


51.0+J 4 


59.7+1.9 ° 




> 48.8 






CI* 


43.2+H 


49.S±| 2 




> 41.1 






Civ 


251+ 27 






> 251 






N v 


109+11 


118+2 5 


142+ 6 3 ; 4 


list 30 - 4 






Ovi 


604±i 




1199+. 272 


sie±2| 5 






Sin 


10 5+ - 6 


13 7+ 12 
'-0.9 


q + 3 . 5 
oo.o_ 2 g 


-,, -.+ 23.1 
-0.9 






Sin* 


i. M +g : i» 


, 1S +0.18 


1 59+ ' 20 
J - aa — Q.17 


1 1»+ - 61 






Sim 








> 9.08 






Si IV 


36.3+.°' 2 


49 4+ 3 7 
^ y -^_2.6 




> 46.8 






Siv 


54.0±*-j 






> 49.5 


T3 


-440 


H i 
Cn 
Cn* 


275t 2 
19. 9+ 2 'o 
7 - 66 tl'.2? 


4230+ 790 b 




4230+ 790 

> 17.9 

> 6.39 






Crv 


301+J 


336+g 


432± 73 


336t 6 69 






N v 


143+j 


167+ 78 




167+ 78 






Ovi 




644+J 7 




644+J 7 






Sin 


< 0.66 






< 0.66 






Sim 


, qo + 0.09 

3.B8_Q Qg 






> 3.80 






Siiv 


13 2+ 11 


16.9±1;S 


ir.3±i;« 


15 9+ 29 
io - — 1.0 






Siv 


< 24.2 






< 24.2 


T4 


-320 


II l 

Crv 


BT.0+S"* 
62.7+ 2 , 


650+iSg b 




550+ 180 
oou -180 

62.7+. 2 , 






Nv 


46.0+ 2 / 2 




67.9+ 4 ; 8 








Ovi 


335±7 


400+ s ° 




400+1° 


T5B 


-130 


II i 

Crv 


399±™ 
795+.* 04 


6010+ 1200 b 

DUJ.U_ 1200 




6010+ 1200 
ouiu_ 1200 

> 784 






Nv 


935± 8 ° 


1035j:" 7 


1469t 284 


1035+. 718 






Ovi 


> 2608_ 22 






> 2586 






Sim 


. _„+o.os 

lre -0.08 






>1.70 






Si IV 


8 49+ 113 


r , ,4-11.7 
1J ' '-1.9 




13 y+11.7 


T5A 


-40 


II l 

Civ 


88-2+n'l 
— 0.5 

97.4± 2 ; 4 


6010 +1200 d 
oulu -1200 




6 [87.7, 7210] 
>95.1 






Nv 


34.7+^ 


36.7tij 




36.7+ 4 j 






Ovi 


> 109_ 4 






> 105 






Sim 


, ,- + 0.06 






>2.19 






Siiv 


5 55+ ' 72 


,„ 7 +176.8 

ij. r_ 2 6 




>13.7 



a) The integrated column densities for the three absorber models. The quoted error arise from photon statistics only and arc computed using the 

technique outlined in I Gab el et al] |2005bh. 

b) Estimates from lpunn et all J 20101 ) . c) Using the covering solution of Siiv(see text) . djT^Hirm et all I2OI0I) do not make the distinction between 
the two sub-components in trough T5. so we report an identical PC value in the shallower trough T5A to be considered as a conservative upper 

limit since the bulk of the column density is coming from T5B. 
e) The lower is fixed by the detection of LylO associated with that component (see text) and the upper limit is given by the absence of a H I 
bound free edge, f) Adopted values for the photoionization study (sec text). 
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and are only lower limits given the saturation of the line profile. Finally, the absence of 
a bound-free edge for H I in the FUSE data allows us to place an upper limit on the H I 
column density of 10 16 ' 2 cm -2 . 

3.2.2. Civ 

Absorption troughs associated with C IV are found in all components of the outflow. 
The velocity range of the outflow (~ 700 km s -1 ) being greater than the separation 
between the components of the doublet Civ AA 1548.200,1550.770 (~ 500 km s" 1 ) limits 
the possibility of deriving a partial covering and power law solution in several components 
of the outflow because of self-blending between the red and blue lines of the doublet. While 
the blue component of trough Tl is free of known blending, its red component is blended 
by the blue component of trough T4. However, we note that the non-blended part of 
trough T4 exhibit the 1:2 strength ratio between the doublet components as expected in the 
AOD model. Assuming that the covering in trough T4 is not a strong function of velocity, 
blending by the red Tl line is limited, suggesting that the trough Tl is also close to AOD. 
Both red sub-components of trough T2 are severely blended by the blue components of 
troughs T5A and T5B. A lower limit on the column density can be placed on trough 
T2 by computing the AOD solution on the non-blended blue line, while a lower limit on 
components T5A and T5-B is placed by the AOD solution on the non-blended red line. 
Kinematic component T3 is the only component where neither C IV of the lines are affected 
by self-blending allowing us to determine the ionic column density using all three absorber 
models. 

3.2.3. Nv 

Absorption troughs associated with Nv are observed in every kinematic component 
of the outflow. Excepting the lower velocity section of trough T4, a high covering fraction 
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Fig. 4. — Comparison between the higher S/N FUSE 2002 and the deconvolved COS 2010 
IRAS F22456-5125 spectra in the rest frame O VI region. While the continuum level in 2010 
is 17% higher than in the 2002 FUSE observation, the absorption trough profile in the red 
component do not show significant changes given the limited S/N of the FUSE observation 
and the larger aperture used by the latter. Changes in the blue line profile seems to be 
observed in the blue component of the Ovi doublet, however, the line profile is located 
at their edge of the COS detector in a region where the S/N is lower and the wavelength 
solution is inaccurate. 
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indicated by the similarities in the column density computed with the AOD and PC 
methods. The higher discrepancy observed for component Tl comes from the fact that 
the blue Nv residual intensity is significantly below the expected level assuming the AOD 
scenario. Since no blend is known to affect the blue Nv, this difference is probably due to 
a slight overestimation of the emission model in that region, so that the column density 
determined on the red N v line is probably a reliable measurement of the ionic column 
density in this component. 



3.2.4. Ovi 

The O VI troughs are located at the edge of the COS detector, where the poor 
wavelength solution and lower S/N limit the constraints we can put on the ionic column 
density. While a higher order correction of the wavelength solution is probably needed at 
the edge of the detector, we use here a single velocity shift of both the red and blue O VI 
lines (respectively 12 and 17 km s _1 ) in order to align the core of the strongest kinematic 
components with the centroids of the Nv A 1238.820 ones. This first order correction 
seems to be sufficient for several components, however the match is not totally convincing 
for others where a clear shift between the centroid of the blue and red line persists (see 
Figure E]). 

Several components of both red and blue O VI lines are blended by known ISM lines, 
further limiting the accuracy of the column density estimates derived for this ion in several 
components. The blue Ovi line is blended in component Tl by a weak Fell A1133.880 
and by Nl A1134.420 and A1134.170 in trough T2 and by Nl A1134.980 in trough T3. 
The troughs of stronger ISM lines associated with Fe II and N I at longer wavelengths are 
shallow, indicating the blends only marginally affect the line profiles of O VI. In trough T2, 
the red line of O VI is affected by an unidentified blend. 

The integrated column densities derived for component Tl using the AOD and PC 
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model reveal a 60% departure from AOD suggesting the PC gives a more realistic estimate 
of the column. The same behavior is observed in component T2. The apparent optical 
depth ratio between the red and blue lines in components T3 and T4 is close to the expected 
laboratory value, suggesting, like in the case of N v, that the AOD determinations can be 
used. Trough T5B is totally saturated in its core and significantly saturated in the wings 
(I r ~ ib). The AOD determination in this case is a strict lower limit on the column density 
since even parts of the red line profile present a residual intensity close to zero (r ^> 1). The 
finite ionic column reported in Table [2] for component T5B is derived using a maximum 
optical depth of r = 4 for these velocity bins. Finally trough T5A shows some evidence 
for partial covering effects. However, the presence of an emission hump around km s _1 
suggests a possible underestimation of the emission model at low velocities around the blue 
O VI line, decreasing the apparent departure from the AOD scenario for this component. In 
order to account for this effect we use the AOD measurement on the red line as a lower 
limit on the column density in component T5A. 



3.2.5. Sim 

Si ill signatures are identified in four of the kinematic components (T2, T3, T5A 
and T5B) of the UV outflow. A weak feature associated with component Tl may also 
be detected in the continuum noise (at less than the 2 a level), though, due to the limits 
of deconvolution, the deepness of the feature is close to other ripples observed in the 
continuum and is probably a false detection. Given the nature of the detection, we report 
an upper limit on the ionic column of Si III for this kinematic component. For the other 
kinematic components, we are only able to place a lower limit on the Si ill column given the 
impossibility of deconvolving the optical depth from the covering fraction for singlet lines. 
We however note that the Si ill trough associated with the kinematic component T2 shares 
a residual intensity identical to the one observed Si IV blue (see Figures [3] and [5]) . Given 
the non-black saturation noted in Section 13.2.61 in component T2 for the Si IV line, this 
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Fig. 5. — The observation of absorption troughs associated with three different ionization 
stages of the Silicon atom allows us to better constrain the photo-ionization analysis of 
the absorber. We also note that the residual intensity in the core of the Si ill T2 troughs 
match the saturated blue Siiv suggesting a non-black saturation of the Si ill (see text). This 
behavior is also observed in trough T3, but leading to a different conclusion due to the non 
saturation of Siiv in that case (see text). 
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observation suggests a net saturation of the Si ill trough whose deepness is mainly reflecting 
a partial covering effect. We note that the residual intensity of the non saturated blue Nv 
line is similar to the one observed in the saturated blue Siiv(see Figure E]), which can be 
a coincidence or point to the fact that the PC model does not constitute an appropriate 
model of the absorbing material distribution. The residual flux in component T3 of Si ill 
also shares a similar depth with the blue Si IV line, however in this case the high covering 
fraction deduced from the residual fluxes observed in the Si IV doublet lines prevent us from 
drawing a similar conclusion. 

3.2.6. Si iv 

Si IV troughs are identified in components T2, T3, T5A and T5B. Trough T2 
shows a significant departure from the AOD model suggesting a better description by an 
inhomogeneous model. Such effect is even stronger in component T5A where the red and 
blue Si IV line profiles matches almost perfectly over the whole component, only allowing us 
to derive a conservative lower limit on the column by assuming an optical depth limit of 
t = 4 in the saturated part of the system. Trough T3 shows a high covering as revealed 
by the small difference between the ionic column densities derived by the AOD and PC 
method. A similar behavior is also observed in trough T5B but with a higher discrepancy in 
the columns due to the difficulty of getting reliable PC measurements for shallow troughs. 

3.2.7. SlV 

S IV is observed in kinematic component T2 as a shallow trough. The ionic column 
density is estimated by applying the AOD method. A shallow feature (less than the 1 a 
level) in the normalized continuum coincides with the expected position of the S IV line in 
trough T3 (see Figure [3]), but the feature is similar in depth to other patterns observed in 
regions free of lines. For this reason we report an upper limit on the S IV line in trough T3 
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by fitting the template of the blue Si IV line to the la noise in that region. 



3.2.8. The density diagnostic lines 

Absorption troughs associated with excited states of Sill and Cll are observed in 
kinematic components T2 (both Si II and C il) and T3 (only C il) allowing us to determine 
the number density and hence the distance to the outflowing material at their origin (see 
Section [5]) . 

Four absorption troughs from the Sill resonance line (A 1190.42, 1193.28, 1260.37, 
1526.72) free of obvious contamination are identified within the COS range associated with 
the kinematic component T2 of the outflow. The weaker A 1304.37 transition (detected 
at less than the 2 a level) is located in a region of lower S/N|j, and is barely detected at 
the S/N level in that region. Observation of four lines emanating from the same state and 
spanning a range of oscillator strengths allows us to further investigate the absorber model 
by over-constraining the residual intensity equations. In this case, we have four residual 
intensities to be fit by two parameters (in the PC and PL models). However, evaluation of 
the fits to the data by these models requires the knowledge of reliable oscillator strengths 
of each line. 

Despite a number of theoretical studies, l arge uncertainties re main in the computed 



Bautista et al 



2009 



for details). Using the 



oscillator strengths of the Si II transitions (see 
oscillator strengths from NIST for the quoted transitions (rated either B+ or C in the 
database), we find that the relative strength order of the lines matches the observed residual 
flux for the A 1190.42, 1193.28 and 1260.37 lines and the weak detection of the A 1304.37 
transition. This is not the case for the A 1526.72 line, which is expected to be weaker than 



3 Given the redshift of the IRAS F22456-5125, the Sill A 1304.37 transition is located in 
a spectral region at the red edge of the G130M grating range and at the blue edge of the 
G160M grating range. 



-23 - 



the A 1190.42 and A 1193.28 lines (hence having a larger residual flux), but for which we 
observe a smaller residual intensity across t he trough. The proble m persists when using 



the updated oscillator strengths reported in 



Bautista et al. 



fl2009f ). While this could be 



due to a blend, the narrowness of the trough and its location away from any known ISM 
lines does not support this scenario. For this reason, we use only the A 1190.42, 1193.28 
and 1260.37 resonance lines to compute the column density. We present in Figure [6] results 
of the simultaneous fits of the three Si II lines performed using the PC and PL absorber 
models. The PC model reveals a small covering factor across the trough, ranging from 0.4 
in the wings to 0.5 in the core. A better fit to the observed line profiles is provided using 
the power law model. The derived power law exponent is close to a = 10 across the trough, 
corresponding to a peaked distribution suggesting that more than half of the emission 
source is actually covered by optically thin material with r < 0.1. The results hold if we 
introduce the weak A 1304.37 transition in the computation, leading to changes in column 
densities that are less than 10% for both PC and PL absorber model. The column density 
derived using the PL absorber model is 2.5 times larger than the one assuming the PC 
model, potentially suggesting an underestimation of the column density when using the 
PC model. However, the PL method can overestimate the column density since a good 
fraction of it originates in the optically thick regi on (V > 4) of the distri bution, to which 



the residual flux is not sensitive (see discussion in 



Arav et al. 



2005 



20081). Finally we note 



that the integrated column density over T2 only diff ers by less than 7 % w hen one uses the 



oscillator strength from NIST or the one reported in 



Bautista et al 



( 120091 ) 



The strongest transition (A 1264.73) associated with the excited Sill (E=287 cm -1 , 
Sill* hereafter) is firmly identified in kinematic component T2 (see Figure |5]). A 
shallower absorption feature corresponding to the weaker excited transition (A 1194.50) is 
distinguished at the S/N of the COS observations. Detection of two transitions from an 
ion with the same low energy level allows us, in principle, to derive a velocity dependent 
solution of the column density using the PC or PL method. However the shallowness of the 
troughs coupled to the limited signal to noise prevent us from computing a reliable solution 
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Fig. 6. — Partial Covering (left) and power law (right) model fit of three Sill resonance 
lines associated with the kinematic component T2 of the UV outflow of IRAS F22456-5125. 



The original data and their error are plo tted in solid whi 



dotted line. The reduced x 2 value (see 



Arav et al. 



e the fitted fluxes are plotted in 



20081 ) of the simultaneous fit for the 



three transitions is given in the bottom right of each panel. The Covering solution C(v) 
is only presented in regions where the residual flux in the weakest lines does not affect the 
determination. 
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across the trough. Nevertheless we can still derive the column density associated with Si II* 
assuming a PC or PL model by using the velocity dependent solution of the covering factor 
C(v) or the power law parameter a(v) derived above from the fitting of the resonance level 
transitions of the same ion. With either C(v) or a(v) fixed in the equations of the residual 
intensity, we constrain the set of 2 equations for the observed residual intensities in the 
Si II* lines and are able to derive the velocity dependent column density solution for both 
absorber models. In trough T3, we put an upper limit on iV(Siii) < 0.66 10 12 cm -2 due to 
non-detection (less that the 1 a level) of the stronger A1260.37 and A1264.73 lines in the 
COS spectrum. 

Cn A 1334.53 (E=0 cm" 1 ) and the blend of Cn A 1335.66,1335.71 (E=63 cm" 1 , 
hereafter ClI*) are detected in components T2 and T3 of the outflow. Having only one 
line for each lower level does not allow us to deconvolve the effects of partial covering and 
population ratio of the level, allowing us to only provide an AOD estimate. In the case of 
Si II we saw that the covering derived for that ion was quite small, having a covering close 
to 0.5 in the core of the trough. Looking at the residual intensity in the core of the C II line 
profile in Figure [TJ it is clear that the covering of that line is larger than 0.5 Thus using the 
Si II covering solution does not allow us to reproduce the observed C II profiles. In order to 
tentatively estimate the effect of covering on the C II and C II* columns, we compute the 
ionic column density using the covering solution derived from Si IV, a medium ionization 
species, and report it in Table |2j While we observe a small increase of the derived columns 
using this PC model, the ratio of column density between the resonance and excited states 
remains identical (as expected given the similar residual flux inside the C II and C II* 
troughs), strengthening the density diagnostic obtained from these lines. 

4. Photoionization analysis of the absorbers 

In order to derive the physical properties of each kinematic component of the outflow, 
we solve the ionization equilibrium equations using version c08.00 of the spectral synthesis 
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Fig. 7. — Line profile of C II and C II* rebinned to a common 5 km s 1 resolution velocity 
scale. Absorption troughs are identified in kinematic components T2 and T3 of the outflow. 
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code Cloudy (last described by iFerland et al.lll998l ). We model each absorber by a 
plane-parallel slab of gas of constant hydrogen number density (n H ) and assume solar 
elemental abund ances as given in C loudy. The spectral energy distribution (SED) we use 



was described in 



Dunn et al. 



(l2010l ). Using the grid-model approach described in Paper I, 
we find a combination of total hydrogen column density (Nh) and ionization parameter that 
best reproduces the observed ionic column densities reported in Table [21 The ionization 
parameter depends on the distance (R) to the absorber from the central source and is given 
by 

Qh 



Uh = r ^- > (3) 

where Qh = 2.5 x 10 55 s _1 is the rate of hydrogen ionizing photons emitted by the object, 
and c is the speed of light. We estimate the hydrogen ionizing rate Qh (and also the 
bolometric luminosity Lb i) by matching the flux of the model SED to the de-reddened 
observed flux at 1150 A(rest-frame) using a standard cosmology (H =73.0 km s _1 Mpc -1 , 
tt A =0.73, tt m =0.27). 



The COS observations sho w a wealth o : 
FUSE observations discussed in Dunn et al. 



absorption lines compared to the earlier 



(l2010l ). This allows us to derive more accurate 



physical properties of the absorbing clouds associated with the UV outflow. In the following 
subsections, we describe the photoionization solution derived for each kinematic component. 
As discussed in Section 13.2.1} the physical properties of the absorber do not appear to 
change between the FUSE and COS e pochs. Therefore, we use the column densities of Hi 



and C ill derived in iDunn et al. 



(120 10h from FUSE data. 



We characterize the maximum likelihood for the model of each kinematic component 
by the merit function : 

2 



X 



E 



log N i<mod - log N { 



obs 



(4) 



log N itQbs - log (N i>obs ± a, 

where, for ion i, N ij0bs and N^ mod are the observed and modeled column densities, 
respectively, and <jj is the error in the measured column density. We prefer this formalism 
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to the traditional definition of x 2 since it preserves the multiplicative nature of the errors 
when dealing with logarithmic values. 

4.1. Troughs T2 and T3 

The physical parameters of component T2 are constrained by ten ionic column densities, 
eight from COS data along with Hi and Cm from FUSE data (keeping in mind that the 
latter have been obtained at a different epoch). The ions span a wide range of ionization 
stages from C II and Si II to N v and O VI. A plot of the results for a grid of photoionization 
models for trough T2 is presented in Figure |HJ The least-squares single-ionization parameter 
solution (x 2 = 1147) is marked with a square in the Njj — Un plane, and predicted values 
for ionic column densities are given in Table [3j The C II and Si II column densities predicted 
by that model are underestimating the observed column densities by one and two order of 
magnitude, respectively, therefore making the model physically unacceptable. Due to the 
poor fit of this model to the data, we use a two- ionization component model (x 2 = 6.7), 
which is depicted by diamonds in Figure [HJ All of the ions from the COS data are fit well 
with the two-component model (see Table EJ). 

For component T3, we have column density measurements for seven ions in the COS 
spectrum, along with Hi and Cm from FUSE data and an upper limit on Sill due to 
non-detection of the stronger lines in the COS spectrum (see Section f3.2.8p . A grid- model 
for trough T3 is plotted in Figure [9j A single-ionization parameter solution (x 2 = 1008) is 
marked with a square in the Nh — Uh plane. This solution fits all the lines within a factor 
of ~ 3 (see Table Hj). An improvement to the fit (x 2 = 5.21) for most ions is provided by the 
two-ionization parameter solution marked with diamonds in Figure with the exception 
being an over-prediction of Sill by a factor £ 4 (see Table HJ). 
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Fig. 8. — Photoionization modeling of kinematic component T2. The plotted lines repre- 
sent slab models whose predicted Ni on matches the estimated values. Solid lines indicate 
a measured column density, dotted dashed lines a lower limit on the column density and 
dotted lines an upper limit on the column density. The error bars due to the photon noise as 
well as the systematic uncertainties in the absorber model are represented as a shaded area 
for each ions when an estimate is available. The black diamonds marks the two ionization 
component model that best fits the estimated N ion while the black square marks our best 

«incrlp inni^ntinn mnrlpl 
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Table 3. Photoionization models for component T2 



Ion log(N ton ) (cm- 2 ) log(f^f) log N mod log(f^f) 
Adopted a SI b Tib TI b. Tl lo +Tl hl 



log f/ff 




-1.9 


-2.8 


-1.3 




log 




19.1 


18.8 


18.7 


19.1 


Hi 


G [ 15.97, 16.19 ] 


-0.25 


16.38 


14.66 


+0.20 


Cii 


> 13.95 


-1.16 


14.13 


10.76 


+0.18 


Civ 


> 14.40 


+0.54 


14.32 


13.96 


+0.08 


Nv 


14 07+ 10 


+0.40 


12.80 


14.05 


0.00 


Ovi 


14 01 +0-25 


-0.18 


12.18 


14.91 


0.00 


Sin 


10 i 7 +0.37 
io - 1 ' -0.03 


-1.91 


13.23 


8.45 


+0.06 


Sim 


> 12.96 


-0.18 


14.01 


10.25 


+1.05 


Si iv 


> 13.67 


-0.37 


13.85 


11.14 


+0.18 


Siv 


> 13.69 


+0.11 


13.59 


12.05 


-0.09 



a) Adopted column densities reported in Table [5] b) The label SI corresponds to the single ionization model 
while TI; and TIftj are the low and high ionization phase of the two ionization model of the absorber. 
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Table 4. Photoionization models for component T3 



Ion log(AT ion ) (cm- 2 ) bg(^f ) log N mod log(f^f) 
Adopted a SI b Tib TI b. TI io +TI w 



log U H ■■■ -2.0 -2.7 -1.7 

logA^ ••• 19.0 18.2 18.8 18.90 

Hi 15.63t^g + - 10 15 - 69 15 - 20 +°- 19 

Cn > 13.39 -0.45 13.42 11.97 +0.05 

Civ U.53t°oH +0.36 13.85 14.49 +0.05 

Nv U.22±° 2 l +0.08 12.44 14.25 +0.04 

Ovi 14.81±g;^ -0.35 11.93 14.72 -0.09 

Sin <11.80 c -0.29 12.47 10.23 +0.67 

Sim > 12.58 +0.38 13.36 11.84 +0.79 

Si iv 13.67t°;^ +0.21 13.27 12.49 +0.13 

Siv < 13.38 +0.44 13.06 13.18 -0.04 



a) Adopted column densities reported in Table [2] b) The label SI corresponds to the single ionization 
model while TI/ and TI^ are the low and high ionization phase of the two ionization model of the 
absorber, c) Upper limit set by the non detection of the stronger Sin and Sin* in that component. 
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Fig. 9. — Photoionization solutions to trough T3. Similar presentation as Figure EJ 
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Fig. 10. — Photoionization modeling of kinematic component Tl. Identical description to 
Figure El 




Fig. 11. — Photoionization modeling of kinematic component T4 
Figure El 



Identical description to 
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Fig. 12. — Photoionization modeling of kinematic component T5B 
Figure El 



Identical description to 
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Fig. 13. — Photoionization modeling of kinematic component T5A Identical description to 
Figure [HI 
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4.2. Troughs Tl, T4, T5A and T5E 

For kinematic component Tl, we essentially have five constraints (ClV, Nv, Ovi, 
Sim and Hi) defining the region of the (Nh,Uh) parameter space able to reproduce the 
estimated ionic column densities. Visual inspection of Figure [TU] suggests a solu tion around 
l og Uh ~ —0.8 and log Nh ~ 18.6, consistent with the upper limits on Cm from 



Dunn et al. 



(120101 ). This solution (x 2 = 86) accounts for every constraint to within 0.25 dex. A 
better solution can be found by relaxing the constraint of solar metallicity. Considering 
the scaling of elemental abundances of C, N and O as a function of the metallicity 



Z (IHamann fc Ferland 



1993 



Hamann 



19971 ). an improved solution (x = 0.5) is found for 
a gas of sub-solar metallicity ([Z/Z Q ] ~ -0.4) using an identical log Uh and total hydrogen 
column of log Nh ~ 19.1 for the slab. 

The constraints on the (Nh,Uh) parameter space for trough T4 are presented in 
Figure [TTJ A solution consistent with the measured ionic column densities is found near 
log Uh — —1.4 and log Nh — 18.4 (x 2 = 2.9). The solution is suggesting roughly solar 
abundances of the gas though the slight discrepancy between the measured H I column 
density and that predicted by the solution may suggest a sub-solar metallicity medium. 

The photoionization solution derived for trough T5B and T5A are presented in 
Figures [12] and [13], respectively. Inspection of Figure [12] suggests a solution around 
log Uh ~ —1.6 and log Nh ~ 19.5 (x 2 = 3.6) for component T5B. A least-square fit to the 
constraints available for trough T5A (Figure [T3|) provides a solution near logUn — —2.0 
and log Nh — 18.6 (x 2 = 27). While the saturation observed in the troughs of several ions 
limits the analysis of the physical properties of the gas, the estimated (Nh,Uh) solution is 
able to reproduce most of the ionic columns to within a factor of 2. 
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Table 5. Ionization solution for each kinematic component 



Component 



log([/ ff ) 8 



\og(N H ) a (cm" 





This work 


DvmnlO b 


This work 


DunnlO b 


Tl 


-0 8 +0 - 24 

u -°-0.20 


-1.15 


18.6185 


18.6 


T2 SI C 


w -0.15 


-1.30 


1Q 1+0-26 
0.24 


19.7 


T2 Tlf 


O+0.04 




1 8 S+0.03 
lS.O_0.04 




T2 IT* 


1 O + 0.10 

J-^-O.OQ 




1 8 7+0-07 
i8 -'-0.04 




T3 SF 


9 n+0.15 
_z - u -0.14 


-1.27 


19 0+ 018 

iJ - u -0.18 


19.7 


T3 Tlf c 


9 7+O.O8 
z - ' -0.06 




i o 9+0.07 




T3 Til 


-i 7+O.O3 
x - ' -0.02 




i o O+0.02 
i8 - 8 -0.03 




T4 


, 4 +o.oi 


-1.54 


18 /L+0- 01 

i8 - 4 -o.oi 


18.5 


T5B 


-i £+0.04 
_1 -°-0.04 


-1.05 


io (5+0-06 


20.1 


T5A 


9 n+°-04 
_z - u -0.04 




1 8 R+0.09 

i8 - b -o.io 





a) The error we report on our determinations of Uh and Nh are estimated by size of the contour in the 
(Uh, Njj ) plane of the solutions that have a \ 2 value twice the x 2 of the best fit solution, b) From 



Dunn et al 



(|2010i ). c) The label SI corresponds to the single ionization model, d) TI; Q and TI^ are the 
low and high ionization phase of the two ionization model of the absorber. 
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5. Absorber distance and energetics 

Detection of resonance and excited state transitions from Sill and Cll in troughs 
T2 and T3 allows us to determine the distance to these two kinematic components from 
the central source. As can be seen from the definition of the ionization parameter Uh 
(Equation [3]), knowledge of the hydrogen number density nn for a given Uh and Qh allows 
us to derive the distance R. When an excited state is populated by collisional excitation, 
the population of t hat state compared to the re sonance level depends on the electron 



number density n e ( lOsterbrock fc Ferland 



2006), which is ~ 1.2 nn in highly ionized 



plasma. Note that photoexcitation could also populate the metastable levels, however, with 
an IR flux ~ 0.1 Jy, population of the metastable levels of Cll and Sill is negligible in 
IRAS F22456-5125. 

In trough T2 we observe resonance and excited states from C II and Si II for which 
column densities have been derived in Section 13.2.81 and reported in Table [2] In Figure HH 
we compare computed collisional excitation models for C II and Si II to the measured ratio 
of the column density between the excited and ground state of these two ions. For the 
Si ll*/S II ratio, the large uncertainty comes from using PC and PL measurements of the 
column density associated with Sill* and S II. The Cii*/Cii ratio is consistent with the 
Sill*/Sll ratio. Given the similar but still significantly different residual flux observed in 
the C II* and C II line profile, while the value of the column density associated with C II* 
and C II can change with the different absorber model, their ratio is less affected since both 
column density will scale with a similar factor (see Section T3.2.8p . For this reason we use 
the Cll*/Cll ratio to derive the electron number density of logn e ~ 1.70l j5 for the low 
ionization phase of component T2. Using the derived ionization parameter of that phase, 
this density implies a distance of R ~ 10.3lf 6 kpc where the errors are conservatively 
computed from the n e range allowed by the Sill*/Sill ratio and the error on the ionization 
parameter. We note that the density derived for this component is consistent with the non 
detection of S IV*, only expected to be observed at higher densities. 
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In trough T3, the only excited state we observe is associated with the C II* transitions. 
Comparing the computed collisional excitation models for C II to the measured ratio of the 
column density between the excited and ground states of this ion (see Figure fTS]) we find 
logn e ~ 1.20loio f° r the low ionization gas phase producing the Cll and Cn* troughs. 
Given the photoionization solution for that phase quoted in Section 14.11 the derived electron 
number density imply a distance of R ~ 16.3lfg kpc where once again the errors on R 
reflect the uncertainty in n e and the ionization parameter. Using the single ionization 
solution for component T3 (see Table [5]) leads to a distance estimate reduced by a factor of 
~2 to R ~ 7.3^13 kpc. Note that in the case of T3, we consider the AOD determination 
of the column densities of C II* and C II given the absence of multiple lines allowing to test 
for the absorber model. Similarly to what is observed for the Si II transitions in kinematic 
component T2, the net difference in residual flux between the excited and resonance troughs 
could lead an inhomogeneous absorber model to predict a smaller ratio (a factor of two in 
the case of the Sill*/Sill ratio of trough T2). The Cii*/Cii ratio derived here could then 
be viewed as an upper limit on the true ratio, the latter being possibly overestimated by up 
to a factor of ~ 2, leading to an underestimation of the distance R by a factor of ~ \/2. 

The reported distances are very large compared to the size of the emission regions in 
AGNs, where we e stimate the broad line emission region to be roughly 0.03 pc in scale 



( IKaspi et al. 



20051 ). Similar distances to narrow absorption line system intrinsic to quasars 



and Seyferts exhibiting lines from e xcited states o 
been reported in the literature (e.g. lHamann et al- 



low ionization species have a 



2001 



Hutsemekers et al. 



ready 



20041 Paper 



I). The kinematic structure and timescale deduced from the line profiles lends support to 
the interpretation of these narrow absorption line systems as associated with episodes of 
mass ejection rather than a continuous wind currently (in the AGN rest-frame) emanating 
from the central regions. In the simplest geometrical picture, we can naively estimate the 
thickness of each shell of ejected material by AR = Nh/uh, where nu is the hydrogen 
number density of the low density, high ionization phase, giving a AR < 2 pc for both 
components T2 and T3. This is assuming that the individual kinematic components can 
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Fig. 14. — Density diagnostic for kinematic component T2. In this figure we plot the 
theoretical ratio of the level population of the first excited states of Cll (E = 63 cm -1 ) 
and of Sill (E = 287 cm -1 ) to the level population of the ground state versus the electron 
number density n e for a temperature of 10000 K (the diagnostic is relatively insensitive to 
temperature for temperatures typical in UV absorber). The ratios derived from Cll and 
Sill (black crosses) imply an electron number density logn e ~ 1.7. The uncertainty on the 
derived n e only accounts for the error on the ratio N ionif /N ion . On the top axis, we report 
the corresponding distance as a function of n e considering the ionization parameter of the 
low ionization phase (log Uh = —2.8). 



-42 - 




Fig. 15. — Density diagnostic for kinematic component T3. We plot the theoretical ratio of 
the level population of the first excited states of Cll (E = 63 cm -1 ) to the level population 
of the ground state versus the electron number density n e for a temperature of 10000 K. The 
ratio derived from Cll and Cll* lines is implying an electron number density logn e ~ 1.2. 
On the top axis, we report the corresponding distance as a function of n e considering the 
ionization parameter of the low ionization phase (log Uh = —2.7). 
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be described as a uniform slab having an internal volume filling factors / of unity. This 
situation is nonphysical since for the inferred temperature of the absorbing gas (T ~ 10 4 
K) the velocity width of the outflow (Av ^ 50 km s _1 )is at least ten times larger than 
the sound speed, and therefore the outflowing material cannot be a sonically connected 
entity. The large Av of this highly supersonic outflow is then probably due to bulk motion 
of the absorbing material. Assuming that the outflowing material is not decelerating, we 
can obtain the dynamical timescale of the shell t s = R/v s ~ 20 Myr, where we choose an 
average outflow speed of v s = 500 km s -1 . Over these 20 Myr taken by the shell to reach 
the distance R, it has been expanding at a speed v exp = FWHM ~ v/10 (see Section |2~T|) . 
so that AR ~ 0.1R. We can use this thickness in order to estimate the actual filling factor 
/ of the shell since AR = Nh / (fnn)- Using the Nh and nn of the high ionization phase 
reported for component T2 and T3 we find / ~ 10 ~ 3 . This number is in ag reement with 



20091 ) based on the 



the low filling factor (/ s < 10 — 0.5) reported by (IBlustin fc Fabian! 
comparison of the observed radio flux and predicted flux at 1.4 GHz, though that study 
was focused on objects possessing an optically thick X-ray absorber. 

Therefore, we assume the geometry of the outflowing material to be in the form of a 
thin (AR < 1/2 R), partially-filled shell, for which we can derive the total mass Mjv in each 
kinematic component i by: 

M Ti = AnR^m p N m , (5) 

where /i = 1.4 is the mean atomic mass per proton, m p is the mass of the proton Q is the 
global covering fraction of the outflow, Nh- is the total hydrogen column density for the 
kinematic component. In the case of a two- ionization component model we simply have 
Nh, = Nh , + Nh- u . . We use Vt = 0.5 since outflow s are detected in about 50% of the 

ft % ,LO Z . fLI ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ 



observed Seyfert galaxies (e.g. ICrenshaw et al.ll2003l ). The average mass flow rate is 



obtained by dividing the total mass of the shell by the dynamical timescale Ri/vi and the 
kinetic luminosity is given by E^. = 1/2 M^vf: 

M Tl = AKR£lpLm p N Hi Vi (6) 
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E ki = 2irR i n(Mi p N Hi yi. (7) 

The advantage of these formulations is that they use of the total hydrogen column 
density iV#. directly derived from the photoionization modeling of the trough and thus are 
independent of the filling factor / of the shell or its AR. We list the computed values of R, 
Mt and Ef. for troughs T2 and T3 in Table [6j 

Note that an instantaneous mass flow rate can be defined independently of the 
dynamical timescale of the outflow by using the physical definition of ina = pAv, 
where p is the mass density of the outflowing material traversing the perpendicular surface 
A with a velocity v. Using the geometry described above this formula simplifies to 
Mt { ins = AirR^QpmpnHifvi, where / is the volume filling factor of the shell. This estimation 
is directly dependent on the filling factor / of the shell (or its radial extent AR), a quantity 
which is not well constrained observationally. We however note that, using the definition 
of the filling factor (/ = Nh/{uhAR)), this instantaneous mass flow rate relates to the 
average mass flow rate defined in Equation [7] by the relation M.T iina = MtJ(AR/R). Since 
AR/R ^ 0.1 (c.f. Section EJ), this means that the average mass flow rates, hence the kinetic 
luminosities, reported in Table [5] are lower limits on the instantaneous mass flow rates. 



6. Discussion and conclusions 

We analyzed the physical properties of the UV outflow of IRAS F22456-5125 based on 
high S/N COS observations. The accurate determination of the column densities associated 
with the multitude of ionic species detected in the COS FUV range allowed us to derive 
the physical parameters (Uh, Nh) of each kinematic components of the outflow. The 
detection of absorption lines associated with excited states of the low ionization species 
C II and Si II in two of the kinematic components allowed us to determine the distance to 
the absorbing material from the central emission source. In the case of component T2, 
the density diagnostic derived from the ratio of Si Ii*/Si II agrees with that derived from 
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C ll*/C II putting the absorbing gas at a distance of ~ 10 kpc. For component T3 only 
the C II diagnostic line is observed suggesting a distance of ~ 16 kpc, though that distance 
could be underestimated by 30% in the case of an inhomogeneous absorber model. 



The photoionization solutions we find differ from those found by iDunn et al.l (120101 ) 
using archival FUSE data(see Tabled]). This is due to the limited number of diagnostics 
available in the FUSE data co mpared with COS data. The total hydrogen column densities 



reported in 



Dunn et al 



(120101 ). derived using only Hland Ovi or Cm, are generally 0.6 
dex higher than the one we find in our analysis. Using the ionization timescale (e.g., Krolik 
and Kriss 1995) along with the assumption that the colum n density of Hi did not change 



Dunn et al. 



(120101 ) estimate a 



between two FUSE observations separated by 21 months, 
lower limit on the distance of ~ 20 kpc to all of the kinematic components of the outflow. 
The distances we find for components T2 and T3 are roughly consistent with this value. 
The small discrepancies may be due to the actual lightc urve over the time period being 



different from simple step-function lightcurve assumed in 



Dunn et al. 



(I2010I) 



Despite the large distance and higher velocity of the outflow compared to the one 
analyzed in Paper I the reported kinetic luminosities E^. in Table [6] are not energetically 
significant for AGN feedback. These scenarios generally require kinetic luminosities to 



be of the order of a few 



Scannapieco &: Oh 



2004 . 



tenths to a few percent of the Eddington lum inosity Lem (e.g. 



Di Matteo et al 



2005 



Hopkins fc 



^lvis 



2008 



20101 ) while in the case of 



), we find E ki ~ 10" 5 L 



Edd- 



IRAS F22456-5125 for which L bol /L Edd ~ 0.16 (IDunn et al.l 
Note that this comparison is probably only a lower limit since it does not take into account 
the fact that the outflow probably decelerated and lost energy through shocks on the way 
from the launching region to the actual location Ri. Moreover, studies of UV outflows in 
Seyfert galaxies typically reveal that an associated warm phase of the outflow, which has 
an ionization parameter substentially large r than the high ioniza tion phase we report here 



(see Table |5]), and usually seen in X-r ays ( 



the kinetic luminosity of the outflow ( iGabel et al. 



Crenshaw et al. 



2005a : 



1999) can carry 70%-99% of 



Arav et al. 



20071). 



Dunn et al. 
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( 12010( 1 analyzed ASCA and XMM-Newton spectra of IRAS-F22456-5125 does not reveal 
any evidence for an X-ray warm absorption edge, however, the limited S/N in these data 
can still allow the presence of a warm phase with significant column density. 

Assuming that the gas is in photoionization equilibrium with the central source, 
we showed that the large number of constraints available for the determination of the 
ionization solution of trough T2 (and in a weaker way in T3) reveals that the absorbing 
material can hardly be described by a slab model characterized by a single Uh and Nh- 
Considering a two ionization solution in which low ionization species are mainly produced 
in a phase whose ionization parameter is ~ 1.5 dex smaller than the phase producing the 
high ionization N v and O VI, we are able to obtain a better match to the measured ionic 
column densities for that component (see Table EJ). If, as inferred from the kinematic 
correspondence, the low and high ionization components are located at the same distance 
from the central source, we can derive the density of the high ionization component to be 
n Hhl = {U Hlo /U Hhi )n Hlo ~ n Hlo /30. 

The observation of several absorption lines corresponding to the Sill transition in 
trough T2 allowed us to test the absorber model by over-constraining the set of fitted 
parameter by 3 residual flux measurements. For that transition we find that the PL 



absorber model describes the observed residua' 
similar to what was found by 



Arav et al 



intensities better than the PC model, 



(120081 ) with the modeling of five Fe II troughs in 



the spectrum of QSO 2359-1241. Either way, the fits suggest that only a small fraction of 
the emission source is covered with optically thick material from that low ionization line. In 
the same kinematic component we observe a high covering fraction for the high ionization 
transition (Civ, Nv, Ovi), almost consistent to a full covering of the emission source 
(continuum+BEL+NEL) by those species. Intermediate ionization species like C II, Si ill or 
even Siiv show clear signs of intermediate covering. These observations suggests a model 
where the low ionization phase is formed by relatively small, discrete clumps of denser 
material embedded in a lower density, higher ionization phase as suggested by iHamann 
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( 119981 ) ; iGabel et al.l (j2005al ) . We however note that the two gas phases are not in pressure 
equilibrium questioning the survival of the low ionization clumps in the more homogeneous 
high ionization phase. 

Comparing the properties of the outflow present in IRAS-F22456-5125 and the bona 
fide AGN outflow observed in NGC 3783 reveals a more complex situation. Albeit their 
similar kinetic luminosity, the in-depth study of the absorption lines presents in the UV 
spectrum of NGC 3783 revealed the expected signs of an AGN outflow : line profile 
variability, high velocities (v ~ 1400 km s - 1 ), high densities ( n P ~ 10 4 ) inferring low 



distances (R < 50pc) to the central source (IGabel et al 



2005al ). The large distance, low 



density, low vel ocity material found in IRAS-F22456-5125 is in comparison also typical of 



galactic winds (IVeilleux et al. 



20051 ) . We reported a similar situation in Paper I in the case 
of the outflowing material present in the quasar IRAS F04250-5718. A key question in 
that case is to determine whether the galactic wind is driven by the AGN or by starburst 
activity. While this question has been investigated in the literature, a definite answer is 
often out of reach for AGNs in which the sustain ing conditions for nuclear activity also 



favor starburst activity (e.g. 



Veilleux et al. 



2005, for a review). While we are not able 



to determine whether the material is AGN or starburst driven, the partial covering and 
the densities higher than the one typically observed in the intergalactic medium deduced 
from our analysis suggests that the material is intrinsic to the host galaxy and is hence 
photoionized by the central source. 

We assumed the ionization structure of the outflow is due to radiation from the 
central source. We investigated the possibility that the absorber is collisionally ionized by 
producing grid- models of Nh versus temperature with a fixed ionization parameter of 10~ 5 . 
At temperatures ~ 10 5 K, we reproduce all the metal lines except Nv and Ovi in trough 
T3. By including another, hotter phase (T « 10 5 ' 5 K), all of the metal lines are reproduced. 
However, Hi is underpredicted by high temperature models (a factor of ~100 for ~ 10 4 ' 8 K 
at log Nh = 19). We are therefore lead to the conclusion that the ionization structure of 
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the absorber is dominated by photoionization. 
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Table 6. Physical properties of the two kinematic components T2 and T3. 

log U H log N H \ogn e R M T \ogE k 

cm~ 2 cm -3 kpc M©/yr ergs/s 



rprt 9 o+0.04 i o O+0.03 i 7 n+0.30 

11 -2.»_ 0-05 18.8_ , 04 1-'U_ 0>15 



-1.3j^° 18.7 + °'°I 0.20+ a31 ' a J ' " L6 ' ~ a9 ' 



-0.09 xo - 1 -0.04 w - i,J -0.19 

T-iq 9 7+O.O8 i o O+0.07 i 9n +0.12 

" z -'-0.06 ±o ' z -0.07 ± ' zu -0.10 

1 7 +0.03 1 o O+0.02 n 9n +0.14,a 

_i -'-0.02 ±o -°-0.03 u - zu -0.12 



} 10.3lf3 S- 1 ^ 41.8_ - 2 



} 16.31ft A - l -o 9 7 41.4^ : | 



a - Computed using the low ionization component number density and assuming that the 
high ionization component is located at the same distance from the central source (see 

Section [6]) 
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